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Abstract. We re-explore the Bi- and Multi-Galileon models with trivial asymptotic 
conditions at infinity and show that propagation of superluminal fluctuations is a 
common and unavoidable feature of these theories, unlike previously claimed in the 
literature. We show that all Multi-Galileon theories containing a Cubic Galileon 
term exhibit superluminalities at large distances from a point source, and that even 
if the Cubic Galileon is not present one can always find sensible matter distributions 
in which there are superluminal modes at large distances. In the Bi-Galileon case we 
explicitly show that there are always superluminal modes around a point source even 
if the Cubic Galileon is not present. Finally, we briefly comment on the possibility 
of avoiding superluminalities by modifying the asymptotic conditions at infinity. 
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1 Introduction 



The combination of independent observational sources such as the CMB, supernovae, 
lensing and baryon acoustic oscillations indicate that the Universe is expanding 
at an increasing rate, possibly driven by dark energy, [1]. Despite a significant 
theoretical and observational efforts the physical origin of the accelerated expansion 
is a mystery. Not only its discovery but also the resolution to the puzzle of its origin 
are an unexaggerated merit of Nobel-prizes. Explanatory attempts fall into three 
primary categories. 

The first solution consists of considering a small cosmological constant A with 
a constant energy density giving rise to an effective repulsive force between cosmo- 
logical objects at large distances [2]. If we assume that the cosmological constant 
correspond to the vacuum energy density, then the theoretical expectations for the 
vacuum energy density caused by fluctuating quantum fields, differs from the ob- 
servational bounds on A by up to 120 orders of magnitude. This giant mismatch 
between the theoretically computed high energy density of the vacuum and the low 
observed value remains for more almost a century one of the most challenging puzzles 
in physics and is called the cosmological constant problem. 

The second solution could for instance consist in introducing new dynamical 
degrees of freedom by invoking new fluids T^ u with negative pressure. Quintessence 
is one of the important representatives of this class of modification. The acceleration 
is due to a scalar field whose kinetic energy is small in comparison to its potential 
energy, causing dynamical equation of state with initially negative values [3]. This 
class of theories might as well as the cosmological constant exhibit fine-tuning prob- 
lems. 

Alternatively, the third solution would correspond to explaining the accelera- 
tion of the Universe by changing the geometrical part of Einstein's equations. In 
particular, weakening gravity on cosmological scales might not only be responsible 
for a late-time speed-up of the Hubble expansion but could also tackle the cos- 
mological constant problem. Such scenarios could arise in massive gravity or in 
higher-dimensional frameworks. 

Concerning theories of higher dimensional theories, the Dvali-Gabadadze-Porrati 
(DGP) model has set the stage for large scale modified theories of gravity [4]. In this 
framework our Universe is a three dimensional brane embedded in a five-dimensional 
bulk. In this higher dimensional setup, the effective four-dimensional graviton on 
the brane carries more degrees of freedom, amongst them the scalar degree of free- 
dom carrying most of the physical impact of this changement. This is explained 
by the graviton acquiring a soft mass m which limits its effective range whereas on 
small scales one recovers General Relativity (GR) through the so called Vainshtein 
mechanism: the additional degree of freedom, the helicity-0 mode, is decoupled from 
the gravitational dynamics via nonlinear interactions of the helicity-0 mode itself, 
[5]. This decoupling of the nonlinear helicity-0 mode is manifest in the limit where 
the four and five dimensional Plank scales are sent to infinity and the soft graviton 
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mass m — > while the strong coupling scale A = (Mpim 2 ) 1 / 3 is kept fixed, enabling 
us to treat the usual helicity-2 mode of gravity linearly while the helicity-0 mode 
7r requires non-linear treatment, [6]. The achievement of the DGP model was the 
existence of a self-accelerating solution sourced by the graviton own degree of free- 
dom, the helicity-0 mode. As promising as this was, it was realized soon that this 
branch of solution was plagued by ghost-like instabilities [7-9]. This issue could be 
avoided for instance including Gauss-Bonnet terms in the bulk, [10]. 

Galilean invariant interactions were then introduced to extend the decoupling 
limit of DGP-gravity [11]. This Galileon model relies strongly on the symmetry 
of the helicity-0 mode tt, which is invariant under internal galileon- and shifting 
transformations tt — > tt + + c. This symmetry can be regarded as residuals of 
the 5-dimensional Poincare invariance in induced gravity braneworld models. These 
symmetries, along with the postulate of the absence of ghosts, imposes drastic re- 
striction on the allowed effective Galileon Lagrangians, as there is a total number of 
only five derivative interactions fulfilling these conditions in four dimensions (omit- 
ting the tadpole). In the context of 5-dimensional braneworld models, these Galilean 
invariant derivative interactions appear as consequences of Lovelock invariants [12]. 
For a review on the Galileon see Ref. [13]. The interesting feature of Galileons to 
allow self-accelerating de Sitter solutions while providing a ghost-free theory [11, 14] 
has generated a flurry of investigations, targeted at cosmological models and their 
observational signatures [15], inflationary models [16-19], gravitational lensing [20], 
spherically symmetric solutions in the vicinity of compact sources [21], Binary Pul- 
sars [22, 23], etc.... 

While cosmological models based on Galileon-gravity were commonly restricted 
to spatially flat backgrounds, there are also possibilities of introducing generalisa- 
tions for non-flat models. One approach being a direct covariantization of the de- 
coupling limit [24, 25]. Naive covariantisation, however, can give rise to ghost-like 
terms in the equation of motion, which can be remedied by a unique non-minimal 
coupling between n and the curvature, [26-28]. The explicit formalism was derived 
in Refs. [26, 27] and the resulting covariantized Galileons are also consistent with a 
higher- dimensional setup [12]. Generalizations to maximally symmetric backgrounds 
with a new Galileon symmetry have been constructed successfully in Refs. [29]. 

Another interesting point to mention is that Galileon-type interaction terms 
naturally arise in theories of massive gravity, which has, in addition, been con- 
structed to be ghost-free be it in three dimensions, [30] or for an generalized Fierz- 
Pauli action in four dimensions [31, 32]. 

Finally, another important property of the Galileon interactions is the fact that 
they are not renormalizable. In other words, the Galileon coupling constants can be 
technically natural tuned to any value and remain stable under quantum corrections 
[6, 33-35]. 

Despite the fact that the Galileon exhibits a broad and interesting phenomenol- 
ogy, they witness a potentially worrying phenomenon, namely the fluctuations of 
the Galileon field can propagate superluminally in the regime of interest [11, 36, 37]. 
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Since the Galileon was introduced as an extention of the decoupling limit of DGP- 
gravity, the DGP model shares the same phenomenon of superluminal fluctuations 
[38]. Superluminal modes are not only generic to Galileon [39] but also to massive 
gravity [40-43]. 

Theories with superluminal fluctuations are sick if they also allow for acausality, 
and configurations with Close-Timelike-Curves (CTCs) are present. Nevertheless 
there are cases in which the superluminal fluctuations come with their own metric 
and causal structure, which can be very different to that felt by photons, and the 
causal cones of these fluctuations might even lie outside the causal cones of photons. 
Regardless of all this, the causal structure of the spacetime can be protected [44] 
if there exists one foliation of spacetime into surfaces which can be considered as 
Cauchy surfaces for both metrics. In theories of Galilean invariant interactions it is 
possible to construct CTCs within the naive regime of validity of the effective field 
theory (as is also the case in GR). Nevertheless, as it has been shown in [45], the 
CTCs never arise since the Galileon inevitably becomes infinitely strongly coupled 
implying an infinite amount of backreaction. The backreaction on the background 
for the Galileon field breaks down the effective field theory and forbids the formation 
of the CTC through the backreaction on the spacetime geometry. The setup of 
background solutions with CTCs becomes unstable with an arbitrarily fast decay 
time. As a result, theories of Galilean invariant interactions satisfy a direct analogue 
of Hawking's chronology protection conjecture, see Ref. [45]. 

The single Galileon scalar field theory has been generalized to a multitude 
of interacting Galileon fields whose origin again can be traced back to Lovelock 
invariants in the higher co-dimension bulk, [34], or such as in Cascading Gravity, 
[46, 48, 49]. Furthermore, they have been extended to arbitrary even p-forms whose 
field equations still only contain second derivatives [47]. 

In this paper, we scrutinize the superluminality in Multi-Galileon theories and 
argue that in these models, the existence of the Vainshtein mechanism about a static 
spherically symmetric source comes hand in hand with the existence of superluminal 
modes. The argument goes as follows: 

1. Superluminalities from the Cubic Galileon for a localized point- 
source: We first show that the mere presence of Cubic Galileon interactions 
guaranties the superluminal propagation of modes in either the radial or the 
orthoradial direction far away from a point source. 

This is intrinsic to the fact that for such configurations, at least one field falls 
as 1/r at large distance as expected from the Coulomb potential. For that be- 
haviour, the matrix encoding the temporal perturbations vanishes at infinity 
while the orthoradial and radial perturbations arise with opposite sign. This 
property is independent of the number of Galileon fields present. 

2. Superluminality from the Quartic Galileon for extended static spher- 
ically symmetric sources: Since the previous result is ubiquitous to any 
Cubic Galileon interactions, the only possible way to avoid superluminalities 
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at large distances, is to set all the Cubic Galileon interactions to zero. In that 
case, we show that the Quartic Galileon always lead to some superluminalities 
at large distances in either the radial or the orthoradial direction again when 
considering a sensible extended source. 

3. Superluminality from the Quartic Bi-Galileon for a localized point- 
source: Even if we restrict ourselves to point sources, we show that the Quar- 
tic Galileon always lead to the propagation of at least one superluminal radial 
mode for some range of r. This result contradicts previous claims found in the 
litterature. 

Our result relies crucially on the assumption that 1. the field decay as the 
Coulomb potential at infinity, 2. that no ghost are present and 3. that the 
Vainshtein mechanism is active (i.e. the Quartic Galileon interactions dom- 
inate over quadratic kinetic terms near the source). The derivation of our 
generic result relies on the interplay between the behaviour of the field at 
large and at small distances. 

4. Superluminality from the Cubic Bi-Galileon for a localized point- 
source: Finally, we show that near a localized source superluminalities are 
also present in the radial direction in a theory which includes only the Cubic 
Galileon. 

The rest of the paper is organized as follows: We start in section 2 with a sum- 
mary of the formalism of the Bi- and Multi-Galileon. We then present the analysis 
needed for the study of the propagation of fluctuations around spherical symmetric 
backgrounds in section 3. In section 4 we study the perturbations around the back- 
ground generated by a point mass at large distances from that source. We show 
that there is always one mode which propagates superluminaly whenever the Cubic 
Galileon is present, regardless of the number of Galileons present in the theory. We 
also find that there are sensible source distributions around which there will always 
be a superluminal mode at large distances even if the Cubic Galileon is absent, for 
any number of Galileons. In section 5 we then consider more closely the case of a 
point mass source when the Cubic Galileon is absent. In particular we study the 
short distance behaviour around a point mass background in the Bi-Galileon theory, 
and we find that there is again always a superluminal mode. In the case of vanish- 
ing asymptotic conditions 7r — > the existence of the Vainshtein mechanism comes 
hand in hand with the existence of superluminal modes. This constitutes a No-go 
theorem showing that superluminal modes are generic to Galileon theories. Finally 
in the discussion section, we comment on the only known loophole to this argument, 
which is to have non-vanishing asymptotic conditions for the field at infinity. This 
kind of set up can arise naturally in theories such as massive gravity which reduce 
to the Galileon theory in some limit. 
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In the rest of this work, we adopt the same notation as in [48] and demonstrate 
in this language that superluminalities can never be avoided 1 in a consistent Galileon 
model. 

2 The Bi- and Multi-Galileon model 

In this work we consider the most general Multi-Galileon theory, in four dimensions, 
which is conformally coupled to matter. This model consists of N coupled scalar 
fields, 7Ti, . . . ,7Tat. For simplicity we neglect gravity in our analysis and study the 
theory on Minkowski space-time. Similarly to Galileon theories [11], the Multi- 
Galileon theory is invariant under internal Galilean and shift transformations 7Ti — > 
7Ti + b^x^ + ci, •■■ ,7Tat— j-7Tiv + b^x^ + . Without loss of generality, we couple 
only one of the N scalar fields to the trace of the stress energy tensor 2 , 

C = £ 2 + £ 3 + £i + £ 5 +Tr 1 T, (2.1) 

where the respective Quadratic £ 2 , Cubic £ 3 , Quartic £ 4 and Quintic Galileon £ 5 
interactions are given by 

£„(7Ti, • • • ,7Tjv) = £ mi ,...,m N (2.2) 

TOiH \-m,]y=n— 1 

with 

£mi,...,m N — { a mi,...,m,N^ + ' ' ' + a m 1 ,...,m N 7r N)£m 1 ,...,m N j (2-3) 

where the a™ ™ are the coefficients for the Galileon interactions. Notice that 
this parameterization allow for a lot of redundancy, so not all the mjv are 

meaningful (many of them can be set to zero without loss of generality). Notice 
as well that in this language these coefficients a's are dimensionfull (the dimension 
depends on mi + • • • + mjv). We stick nonetheless to this notation for historical 
reasons, [48]. In this formalism, all the derivative are included in the £ mii .... mAr 
which can be expressed as 

S mi ,..., mN = (mi + • • • + m N )%\ ■ ■ ■ 8^1 (2.4) 
x [(d^m) ■ ■ ■ (d^d^n,)] ■ ■ ■ [(d n ff*ir N ) ■ ■ ■ (d PmN d°^n N )] , 

using the formalism derived in Ref. [27]. 

1 So long as the Galileon interactions dominate near the source (i.e. , as long as the model 
exhibits a viable Vainshtein mechanism) and as long as one considers trivial asymptotic conditions 
at infinity. 

2 One can always rotate the field space 7Ti, • • • , t^n so as to couple only one field to the trace 
of the stress-energy tensor. This would not be possible if more subtle coupling to matter were 
present, for instance 9 M 7r9i / 7rT /iI/ as is the case in Massive Gravity [32]. However this coupling 
cancels at the background level for static spherically symmetric sources. 
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Bi-Galileon: Specializing this to the Bi-Galileon N = 2 is straightforward. The 
analogue to (2.3) for the Bi-Galileon would simply be 



£ 

*— "71 



0<m+n<4 



with 



^ ^ Q"m,n£m,n T and ^ ^ b mn S mn 



0<m+n<4 



0<m+n<4 



(2.5) 



(2.6) 

The equations of motion for the two scalar fields tt\ and 7r 2 are then 



(2.7) 



where the coefficients a m>n and b m ^ n can be expressed in terms of the parameters 
and /3 m , n as 



(m + l)(a m ,n + 



m+l,n—l, 



and 



m— l,n+l 



). (2.8) 



3 Spherical symmetric backgrounds 

In this subsection, we recapitulate the formalism needed to study the superluminality 
of fluctuations about spherical symmetric solutions. For this we split every field into 
a spherically symmetric background configuration 7r°(r) and a fluctuation 5ir(t,f), 



■K n (t, r) = (r) + <5vr n (t, r ) , V n = 1, • • ■ , N , 
and introduce the N- dimensional fluctuation vector in field space, 

5U(t, r)= : 

\5ir N (t,r)J 

At quadratic order in the fluctuations, the Lagrangian can be written as 
A*,...,™ = \d t U .K.dtlL- l -d r Tl . U . d r U - ^d n Il .V.d n U. 



(3.1) 



(3.2) 



(3.3) 



The kinetic matrix K, and the two gradient matrices U and V are defined as follows: 



K = (1 + -<9 r ) (S x + 3S 2 + 6S 3 + 6E 4 ) 
o 

W = Si + 2S 2 + 2S 3 
V = (l + ^ r )W, 



(3.4) 
(3.5) 
(3.6) 
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where the E matrices depend on the spherically symmetric background configuration 
(and are thus functions of r). In this language the n th matrix E n encodes information 
about the (n + l) th order Galileon interactions C n+ i, 



with 



(Q pi . . . Q fa N \ 
u yi J n U VN J n 1 

\ u Vl J n U VN Jn / 



(3.7) 



/n(yiW 



VN{r 



(3.8) 



i=0 



and for each of the Galileon field, we define, 



yn{r) 



«(r) 



(3.9) 



In terms of the matrix U, the background equations of motion are given by 



/ 



r 2 U(it°(r)).d r 



V 



/7r?(r)\ 



7To r 



\ 



J 







(3.10) 



In particular for a point source of mass M = Awm localized at the origin r = 0, we 
have 



(E 1 + 2E 2 (r) + 2E 3 (r)). 



/ yi( r ) \ 





V o y 



(3.11) 



where Ej is independent of y n and is thus simply a constant, E 2 is linear in the y n 
and E 3 is quadratic in the fields. 

Notice that the expressions (3.4, 3.5, 3.6) for the matrices /C, U and V in terms 
of the E n matrices are universal and do not depend on the number N of fields. 

Focus on the Bi-Galileon: 

In the following we restrict our attention to the Bi-Galileon since we will first focus 
on that case and then generalize our results to the Multi-Galileon case. In the 
Bi-Galileon case, the matrices E n are given explicitly as below: 



ft/ill) with £ = 5>uv*" 



(3.12) 
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The functions y and z appearing in the coefficients f° are shortcuts for 
such that the equations of motion for 7r° and 71"° become simply 

777 

n+^f a 2+n) = -, (3.15) 



r 



A fe + 2(/ 2 fe + / 3 b )=0, (3.16) 

where m = M/Att, and M is the mass of the point particle introduced at r = 0. 
More explicitly, in terms of y and z the two equations of motion are given by 

aioy + a iz + 2 (a 20 y 2 + a n yz + a 02 z 2 ) 

TTl 

+ 2 (+a 30 y 3 + a 2 i2/ 2 z + a 12 yz 2 + a 03 2 3 ) = -3, (3.17) 

( ^11 2 r 2^ 

a iy + &01Z + 2 ^— 3/ + 2a 02 i/2; + 6022 J 

+ 2 (a 21 /3y 3 + a 12 y 2 z + 3a 03 yz 2 + & 03 2 3 ) = . (3.18) 

In terms of the parameters and bij, the £1,2,3 matrices can then be expressed 
respectively as 

s - - (:: • ^ 

s = f2a 2 oy + a 11 z a 11 y + 2a 02 z\ (3 20) 

\a n y + 2a 2^ 2a 02 y + 2b 02 z / ' 

s _ f3a 30 y 2 + 2a 2 iyz + a 12 z 2 3a 03 z 2 + 2a 12 zy + a 21 y 2 \ (3 21) 

' \3a 03 z 2 + 2a 12 zy + a 21 y 2 3b 03 z 2 + 6a 03 zy + a l2 y 2 ) ' 

in the Bi-Galileon case. To get these expressions we have used the fact that for m<n 
we have the correspondences £ m ,n = £n,m\^^n - The exclusion of super luminal mode 
propagation implies that the sound speed of both modes along both the radial and 
orthoradial directions be less than or exactly equal to 1. The two sound speeds in 
the radial direction are given by the eigenvalues of the matrix M. r = 1C~ 1 U and the 
two sound speeds along the orthoradial direction are given by the eigenvalues of the 
matrix M.^ = /C _1 V. Therefore the condition for no superluminality is equivalent 
to requiring that all the eigenvalues of both matrices M. T — I and Ain — I be zero or 
negative (and larger than —1), with I the identity matrix. In the following sections 
we study the behavior of the system in two different regimes, in the large and short 
distance regimes and check explicitly there always exists at least one superluminal 
mode in one direction. 
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4 Superluminalities at Large Distances 



Let us start by summarizing the results find in this section: 



• We first show that if at least one Cubic Galileon interaction is present 3 i.e. if 
£3 does not vanish identically, then superluminal propagation is always present 
at large enough distances from a point source. The only way to bypass this 
conclusion is to remove all the Cubic Galileon interactions for all N fields 
£3 = 0, meaning that any ... mjv with n — 1, • • • ,N and m 1 + ■ ■ ■ + m N = 2 
has to vanish exactly (for example in the Bi-Galileon case, this implies 020 = 
On = «02 = ^02 = 0). If the coefficients are merely small, then one can always 
go to large enough distances where the Cubic Galileon dominates over the 
Quartic and Quintic Galileon interactions. 

• Nevertheless even if all the Cubic Galileon terms vanish £3 = 0, we can 
still find perfectly sensible static, spherically symmetric matter distributions 
around which there are superluminalities due to the Quartic Galileon at large 
distances. 

• As a consequence, we will see in this section that as soon as either a Cubic 
or a Quartic Galileon interaction is present one can always construct a sensi- 
ble matter distribution which forces at least one of the iV Galileon fields to 
propagate superluminaly in one direction (either the radial or the orthoradial 
one). 

We note that the Quintic Galileon interactions £5 always vanish at the background 
level around static, spherically symmetric sources, independently of the number of 
fields, so that if one tries to avoid the above conclusions by making both the Cubic 
and the Quartic Galileon vanish, then there is no Vainshtein mechanism at all about 
these configurations. 

4.1 Superluminalities from the Cubic Galileon 

In the Multi-Galileon case, the background equations of motion for a point source 
at r = are given in (3.11). At large distances 4 , this simplify to 



Si 



/ yi( r ) \ 

V2(r) 
\vN(r)J 



( r 3 \ 





\°7 



(4.1) 



3 It can be an interaction involving just one of the N Galileon fields, or an interaction mixing 
different Galileon fields together, the result remains unchanged. 

4 We assume throughout this manuscript trivial asymptotic conditions at infinity which implies 
that the Galileon interactions ought to die out at large distances. The contributions from Si are 
thus the leading ones at large distances. Consistency of the theory requires that det Si 7^ (so 
that the theory does indeed exhibit N degrees of freedom) and the matrix Si is thus invcrtiblc. 



- 10 - 



Recalling that Si is invertible and independent of the field (Si is a constant), this 
implies that to leading order at large distance about a point-source, the fields die 
off as r _1 

y{r) ~ r" 3 + C(r^ 6 ) 7r°(r) ~ r^ 1 + C(r~ 4 ) , (4.2) 

for at least one of the N fields, as expected from the Newtonian inverse square law 
which should be valid at infinity. As a result, at large distances the S matrices 
behave as follows: 

51 = Si (4.3) 

5 2 = ls 2 + O (r- 6 ) (4.4) 

5 3 = ^S 3 + 0(r" 9 ) , (4.5) 

where the 'barred' matrices £1,2,3 are independent of r. 

As a direct result of this scaling, it is trivial to see that at large distances, the 
kinetic and gradient matrices /C, U and V are given by 

/C = Si + + o(l) , (4.6) 

U = Si + ^2 + O (1) (4.7) 
^ = S i-^ 2 + o(l) . (4.8) 

It is apparent that the perturbations at the order -3 in the matrix /C vanish while 
in the matrices U and V they always come with the opposite sign, hence there is 
always a superluminal direction at infinity. These results coincide with what is 
already known in the case of one Galileon [11]. This is intrinsic to the -g behaviour 
at infinity and to the presence of the Cubic Galileon, and is independent of the 
number of fields. 

The only way to bypass this very general result is to require the matrix S 2 
to vanish entirely, which could be for instance achieved by imposing all the Cubic 
Galileon interactions to vanish 5 . In particular, even if some eigenvalues of S 2 vanish, 
the previous result remains unchanged, as long as S 2 has at least one non- vanishing 
eigenvalue which would imply that the associated eigenmode in field space has a 
superluminal direction (either a radial or an orthoradial one). Only if all the eigen- 
values of S 2 vanish can we evade the previous argument, which can be accomplished 



5 At large enough distances the Cubic Galileon would always dominate over the Quartic one 
(assuming trivial asymptotic conditions at infinity), so imposing a hierarchy between the Cubic 
Galileon interactions and the other ones is not sufficient to avoid superluminalitics. All the Cubic 
interactions should be completely absent. 
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by demanding all the coefficients arising from the cubic Galileon interactions to van- 
ish exactly, in other words if the next to leading interactions arise from the Quartic 
Galileon. 

In the very special case where £2 vanishes entirely (i.e. all its eigenvalues are 
identically 0), then the previous argument needs special care. The contribution from 
S 3 implies that (JC~ l U) and (/C -1 V) do not necessarily have opposite sign. In any 
Multi-Galileon theory one can always tune the coefficients of £3 so that the matrix 
E 2 vanishes identically and so the r -3 scaling is not the leading order correction to U 
and V. For example, for the Bi-Galileon if the parameters of the theory are carefully 
chosen so as to satisfy 020 = ^02C 3 , an = 26 2C 2 and a 02 = b 02 c with c = aoi/601 
then U and V vanish identically at 0(r~ 3 ) for a pure point source and the argument 
given above breaks down. However as soon as we consider an extended source with 
energy density going as 1/ r 3_e would revive £2 and the argument would then again 
be the one above. So even for these special coefficients, there is a whole classs of 
otherwise physically sensible solutions which exhibit superluminal propagation. 

In section 5 we consider this case more closely in the Bi-Galileon scenario and 
find that there is still always at least a superluminal mode for some range of r. For 
instance superluminalities unavoidably arise near the origin through the Quartic 
Galileon, unlike what was claimed in [49]. But first we point out that one can very 
easily construct an extended source for which superluminalities are present at large 
distance for the Quartic Galileons just in the same way as they were for the Cubic 
ones. 

4.2 Superluminalities from the Quartic Galileon about an Extended 



When the Cubic Galileon is absent, the presence of superluminalities about a point- 
source is more subtle to prove and will be done explicitly in the next section. Nev- 
ertheless, even if the coefficients in the Cubic Galileon vanish, we can always find a 
background configuration in which we can see superluminalities at large distances 
arise using the same argument as for the Cubic Galileon. In particular we can 
consider a gas of particles with a spatially varying density of the form 



where r characterizes the scale over which the density varies and M controls the 
overall strength of the density profile 6 . 

In this case the asymptotic behavior of the background fields becomes 



Source 




(4.9) 




(4.10) 



6 This matter distribution can always be imagined for some arbitrarily large radius before being 
cut off. 



for all the fields n = 1, • • • , N, and we find once again using Eq. (3.4) that to order 
0(r -3 / 2 ), /C vanishes and U and V have opposite signs, guaranteeing a superluminal 
direction. 

This illustrates the basic reason we expect any theory that exhibits the Vain- 
shtein mechanism to inevitably contain superluminalities when considering trivial 
asymptotic conditions. Every new source configuration gives rise to a new back- 
ground Galileon field configuration. Because the theory must be nonlinear in order 
to have a Vainshtein mechanism, the fluctuations around this background will prop- 
agate on an effective metric determined by the background. Since the sources are 
not constrained by the theory, we are free to choose any source we like, and so we 
have a lot of freedom to change the parameters in this background metric and create 
superluminalit ies . 

5 Quartic Galileon about a point-source 

In the previous section we have shown that superluminalities are ubiquitous in 
Galileon models. No matter the number of field there is no possible choice of param- 
eters that can ever free the theory from superluminal propagation. The argument in 
the previous section was completely generic an independent of the number of fields. 
It only required the behaviour at large distances (when the non-linear Galileon in- 
teractions can be treated perturbatively). 

In what follows we show that even for a point source in the Quartic Bi-Galileon 
model 7 , superluminalities can never be avoided in a consistent model 8 . This result 
is in contradiction with previous results and examples in the literature, but upon 
presentation of this following argument, the previous claims have been reconsidered. 

The philosophy of the argument goes as follows: We analyze the model both at 
large distances in the weak field limit and at short distances from the point source 
where the quartic interactions dominate (as required by the existence of an active 
Vainshtein mechanism). The requirement for stability (in particular the absence of 
ghost) sets some conditions on the parameters of the theory. We then show that 
these conditions are sufficient to imply the presence of superluminal modes near 
the source. We emphasize that this result could not be obtained, should we just 
have focused on the near origin behaviour without knowledge of the field stability 
at infinity. 

5.1 Stability at Large Distances 

To ensure the stability of the fields, the kinetic matrix /C as well as the gradient 
matrices U and V should be positive definite at any point r. At infinity in particular 

7 i.e. we only consider two fields with mixing kinetic terms and Quartic Galileon interactions, 
but no Cubic interactions. 

8 The only requirements are the absence of ghost, the presence of an active Vainshtein mecha- 
nism, and trivial conditions at infinity. 
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these three conditions are equivalent and simply imply that the matrix Si ought to 
be positive definite. In the case of the Bi-Galileon, this implies 

det Ei = a w b i — > and Tr E x = a w + fe i > . (5.1) 

In terms of the coefficients of the quadratic terms, these two conditions imply 

a 2 

aw > and 6 i > — > . (5.2) 

aio 

The behaviour of the fields at large distance from a point source localized at r = 
is determined by the coefficients of the quadratic terms, (or equivalently by Ei), 

y(r) = £ + O (i) , with n = ^-m (5.3) 



*> = ! + o(^). with Z^^-m, (5.4) 



and as expected, we recover a Newtonian inverse square law behavior for each mode 
at infinity, namely d r ir® ~ dj-Tr® ~ t~ 2 . At this stage it is worth to mention that the 
stability condition (5.2) implies that Y\ > 0, which is consistent with the fact that 
the force mediated by the one field 7Ti that coules to matter is attractive. 

5.2 Short Distance Behavior 

We now study the field fluctuations at small distances near the source (i.e. at leading 
order in r, assuming we are well within the Vainshtein region). From the equations 
of motion (3.17), (3.18) after setting the coefficients of the cubic Galileon to zero) 
near the origin, we infer the following expansion 

y(r) = — + y 2 r + 0(r 3 ) (5.5) 
r 

z ( r ) = ^- + Z2r + Q(r 3 ) , (5.6) 
r 

with 

amvl + flfciJ/iZi + a 12 y x z\ + a 03 zf = — (5.7) 

— Vi + a 12 y 1 zi + 3a 03 ?/i2q + Oq 3 z\ = . (5.8) 

Note that the 0(r°) terms vanish since the Cubic Galileon is not present. 
Expanding E 3 in powers of r, we have 

E 3 = Ef + Ef 5 + • • • = is® + + O (r 2 ) , (5.9) 

where the leading order contribution to E 3 is given by 

= / 3a 30 yf + 2a 2 iy\Zi + a 12 zf 3a 03 z 2 + 2a 12 yiZi + a 21 yf\ , . 

3 \3a 03 zf + 2ax2yiZ\ + ct 2 i?/ 2 36 03 2 2 + 6a 03 ?/i2i + a 12 yf J 
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Solving the equation of motion (5.8) for b 03 gives b 03 = (—a2iy 3 —3ai2y 2 z—9a 03 yz 2 )/(3z 3 ). 
Similarly solving the equation of motion (5.7) for a 30 yields 030 = ((m — 8a 2 iTrr 3 y 2 z — 
8ai27ir 3 yz 2 — Sao 3 nr 3 z 3 )/ (8nr 3 y 3 )). Using these expressions for b 03 ,a 30 and intro- 
ducing the combination B defined as follows: 

B = 3a 03 z 2 + 2a 12 y%z 1 + a 2 \y\ , (5.11) 

we can then write S3 simply as : 

ff = fe _L) • < 5 ' 12 > 

5.3 Stability at Short Distances 

As mentioned previously, me must ensure that the eigenvalues of /C are strictly 
positive. At small distances near the source, the matrix K, can be expressed as 

K = + Oil) . (5.13) 

In terms of yi, z\ and B, the absence of ghost near the origin implies the following 
conditions on the parameters 

detE 3 = -5— > (5.14) 
TVE 3 = 2^1^M±5) >0 , (5.15) 



which are equivalent to 



B 

yi>0, and — <0. (5.16) 



We now use the stability conditions derives at both large and small distances to 
deduce the behaviour of the radial sound speed in this model. 

5.4 Sound Speed near the Source 

Similarly as we did at large distances, we can now compute the 'radial sound speed' 
matrix Ai r = Kr x U near the origin, 

M r = I — 2r 2 (S 3 °)- 1 Ef } + C(r 4 ) . (5.17) 

We note that unlike the Cubic Galileon case described in more detail below, the 
leading order behaviour of Ai is not manifestly superluminal. However, this is not 
enough to guarantee the absence of superluminal modes, we must carefully check 
the sign of the small 0(r 2 ) correction term before making any conclusions. A simple 
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formulation for the matrix £ 3 ^ is given in (5.12), and a similar expression for £3"^ 
can be found in an analogous way, 

3 " V 2 C -hi-ao^-^J ' (5 ' 18j 
with the notation: 

C = 0212/12/2 + 0122/2^1 + a i2?/i^2 + 3a 03 ^2 • (5.19) 

This allows us to compute the radial sound speed 9 : 

c 2 s± = l + r 2 (a' ±VV) + 0(r 4 ), (5.20) 

with a' and V some coefficients that depend on y X) z x , B, m and (a, b)^. So for both 
modes to be subluminal along the radial direction, the following conditions should 
be satisfied: 

a' < 0, V > and a' 2 - V > . (5.21) 

However as we shall see, these are not consistent with the stability conditions es- 
tablished previously. 

The explicit form of the coefficients a' and 6' 10 is given by : 



a 1 



12 

a 



-■5^- (3m(aoiJ/i + 2Cyi + b mZl ) - 2i3T?) (5.22) 

8 / \ 
b' = - 3m g y (( a io2/i + aoiZi)(aoiVi + boxZi) + 2(V) , (5.23) 

with the notation 

D = axoVi + 2a i2iyi + b m z\ . (5.24) 
We may re-express the quantity T> as follows 

V = aioVi + 2aoi^iJ/i + b iz 2 (5.25) 

/ a 01 \ 2 z 2 
= aw ( Vi + —Zi J + — (a w b 01 - a 2 m ) > . (5.26) 
V aio / oio 

Recall from the expression of the kinetic matrix K, at infinity, that the following two 
conditions should be satisfied, (5.1): a 10 > and (eiio&oi — a oi) > 0' which directly 



9 We need to work to 0(r 4 ) inside the square root to get this expression, because we need to 
square the 0(r 2 ) corrections we have calculated. One might worry that the calculation we have 
done is not consistent because we have not worked to C(r 4 ), however one can check that the 0(r 4 ) 
corrections we have neglected here cancel identically and do not contribute to c 2 s at the order we 
are interested in. 

10 It easy to check that b' is always positive. 
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implies that T> is strictly positive. Knowing this, we check whether or not a' < 
and a' 2 — b' > which if true, would imply that both modes are sub-luminal. 

We start with the requirement that a' < 0. This implies that 3m(a iyi + 2(y 1 + 
boiZi) — 2B(a w yf + 2a^\Z\y\ + feoi^i) has the same sign as B. Once this condition is 
satisfied, we check the sign of a' 2 — b' > 0. This quantity is positive only if J 7 has 
the opposite sign as B, where 

? = (a-ioVi + a 01 z 1 )(a 01 y 1 + fo i^i) + 2CV . (5.27) 

In what follows, we will start by assuming that z\ is positive and show that in that 
case the condition to avoid any super-luminal modes cannot be satisfied. The same 
remains true if Z\ is assumed to be negative. We can therefore conclude that the 
Quartic Bi-Galileon interactions always produce a superluminal mode already in the 
configuration about a point source. 

We recall from eq. (5.16) that if z\ > 0, the absence of ghost-like modes near 
the origin imposes the condition B < 0. Furthermore, knowing that T> = aioyf + 
2a>oiZiyi + feoi^i > 0, we can infer that a' negative only if 

2BV 

a m y 1 + 2(y l + b m z 1 <- — <0. (5.28) 

3m 

Then using the fact that T> = a^yf + 2a$\Z\y\ + 601^1 > 0, this implies (knowing 
from (5.16) that y x > 0): 

1 z 

ooii/i + 2( Vl + 6 i2i < C < -^( a oi + b 01 — ). (5.29) 

^ Hi 

Finally to avoid any superluminal mode, the quantity a' 2 — b' should also be positive. 
Since in this case B is negative, a' 2 — b' has the same sign as J 7 , where 

F = (aioZ/i + aoi*i)(aoi!/i + boiZi) + 2(1) 

< (aioZ/i + aoiZi)(a iyi + 601 Zi) - (aoi + b i—)V 

yi 

< -— (aoij/i + boiZi) 2 . (5.30) 

yi 

Since y\ > and z\ > this implies that J 7 < 0. Since a' 2 — b' has the same sign 
as J 7 , we can therefore conclude that if we assume Z\ to be positive and a! < 0, the 
quantity (a' 2 — b') will also be negative, or in other words, there is one superluminal 
mode. This argument was made assuming z\ > 0, however it is straightforward to 
reproduce the same argument for negative z\. If we choose for instance negative 
z\ {z\ < 0) then the condition coming from the absence of ghost-like instabilities 
eq. (5.16) will require this time the opposite sign for B, namely B > and therefore 
T will be a positive number T < — ^ (a m yi + b 01 z\) 2 while the expression (a' 2 — b') 
in eq. (5.23) will have the opposite sign to J 7 and therefore again there would not 
be any choice of coefficients (a, b)ij to make both modes (sub)luminal. With this we 
have proven that there is no generic choice for the parameters and b^ near the 
origin which would prevent the propagation of superluminal modes. 
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5.5 Special Case of Dominant First Order Corrections 

In the previous section we proved that close to the source there is always one mode 
which propagates superluminaly in a generic theory where only the Quartic Galileon 
is present. However we made a technical assumption in Eq. (5.17) that Eg was 
invertible, or equivalently that we did not make the special choice of parameters 
ay, boi which gives B = (implying that the leading order pieces in E 3 were strictly 
larger than the first order corrections). However we could consider a specific choice 
for which some of the leading order pieces of £3 vanish and the subleading pieces 
become dominant. Therefore in this section we will examine this possible loophole 
more closely. We will find that even in this case one always finds that a superluminal 
mode is present. When B = 0, E3 takes the following trivial form: 

( 3m q \ 
2 ooj- (5 ' 31) 

The stability of the theory now depends not only on the leading behavior of the 
kinetic /C and radial derivative U matrices, but also on the subleading behavior. In 
terms of the E matrices, /C and U take the following form 

E (0 £ (0 

/C = 2^ r + £ 1 + 6£i ni) , W = 2^- + E 1 + 2E^°. (5.32) 

The theory is stable only if K. and U have positive eigenvalues, or in other words 
only if the following quantities three quantities are positive: 

yi > 0, 

A? = 2 (6(a i + 2Qy lZl + 5&oi*?) > 0, (5.33) 

A' ee 2 ( 2 (°oi + 2C)ViZi + boizf) > 0. 

Now we construct again the radial speed of sound matrix A4 r = K,~ l lA in this specific 
case with B = 0. We can write the trace and determinant as 

trM r = + +r 2 r + 0(r 4 ), 



A? 



detMr = + r 2 5 + C(r 4 ). (5.34) 
Ai 

where r and 5 are functions of the given parameters (however we will only need r — 5 
as shown below). The speed of sound is given by, to G(r 2 ), 

• If > X 2 2 
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In this case we will have superluminal propagation if and only if r — 5 > 0. We show 
that one always has r — 5 > in this case by carefully making use of the stability 
constraints in Appendix A. 

• If A? < \ 2 2 

C±= \l-^-^ (5 ' 36) 

Since A^ < Af, in this case superluminal propagation is guaranteed. 

One might argue that this only guarantees superluminalitiy at the origin which 
is inside the redressed strong coupling radius of the theory, where we can no longer 
trust the results of the theory. However, we note that explicitly factoring out powers 
of M and A that the speed of sound is given by 

^f+G^A^- (5 - 37) 

where ry = (M/Mp{)^ 3 /A is the Vainshtein radius and where f and 5 are dimen- 
sionless. Since ry > A -1 , and since the redressed strong coupling radius is always 
smaller than A -1 , there is a range of r in which we can trust the theory and we can 
also trust the leading order behavior of the speed of sound above. 

6 Cubic Lagrangian near the Source 

Lets have a quick look into the contributions coming from a Cubic Bi-Galileon 
theory 11 near the origin and study the superluminality. In the section 4.1 we had 
seen that the existence of the Cubic Galileon guarantees superluminal propagation 
at infinity. Now lets also see the effect of a pure Cubic Galileon term on short 
distances. We quickly recall the equations of motion in the Cubic Galileon case near 
the origin here again: 



2 (a 02 z 2 + a n yz + a 20 y 2 ) = — (6.1) 



2 

At short distance the fields then behave as 



2 ( b Q2 z 2 + 2a 02 yz + \a liy 2 ) = . (6.2) 



y( r ) = ^ + y^ + °^ 12 ) > ( 6 - 3 ) 



z{r) = ^ + z 2 + 0(r 3 / 2 ) . (6.4) 



ii 



i.e. a Bi-Galileon theory where only the Cubic interactions are present 
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The leading order matrix E^ can be expressed as follows (after use of the equations 
of motion): 



— — —C c 

yi yi 

C -*±C 



(6.5) 



with the notation 



C = a n yi + 2a 02 2i 



(6.6) 



The stability condition in the short distance regime requires det/C ~ detE^ > 
and Tr/C « TrE^ > 0: 



TrE 



(/) 



m — Czi 
Vi 



> 



detE 



(0 



-mC 



> 0. 



(6.7) 
(6.8) 



These conditions imply — < and y\ > 0. After using the next-to-leading order 
equations of motion to simplify the next-to- leading order matrix E^ , we find 



with (3 



a 2Z2 + 



aio „ zi QRZi 



22/i 

2/3 



yi 



&oi 
2 



2/3 

2zi 



2(3 



(6.9) 



Assuming E^ is invertible (which is the case if there is at least one non- 
vanishing Cubic Galieon interaction), the matrix Ai r is given by 



M r 



KT X U 
4 2 



3E 



(0 



.3/2 



2r 3/2 



v (0 

^2 



Si+4E^ 



(nl) 



.3/2 



:Sx + 2E 



(6.10) 

(6.11) 



This in turn implies that the Cubic Galileon also gives rise to superluminal propa- 
gation near the origin. If on the other hand we consider the possible loophole with 
vanishing determinant of the leading matrix E^ (choosing parameters such C — 0) 
nothing changes. The matrix Ai r has still one eigenvalue going as 4/3 + 0(r), 
and another eigenvalue whose leading behavior depends on the signs and relative 
strengths of and Z\. But the existence of one eigenvalue that is 4/3 at leading 
order is enough to prove the existence of superluminalities in that regime as well. 



7 Discussion 

In this paper we have shown that Multi-Galileon theories inevitably contain super- 
luminal modes around some backgrounds, for any number of Galileon fields. At 
large distances from a static point source, we have shown that if the Cubic Galileon 
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is present (even if its coefficients are very small), it will eventually dominate over 
the other Galileons and lead to a superluminal mode. Even if no Cubic Galileon 
interactions are present (i.e. all the Cubic Galileon coefficients are exactly zero), we 
find that there are simple, perfectly valid matter distributions (such as a static gas 
of particles whose density falls of as r -3 / 2 ) around which perturbations propagate 
superluminally. 

We also considered the case studied in the journal version of Ref. [49] of per- 
turbations around a static point source in the Bi-Galileon when only the Quartic 
Galileon is present. By studying the speed of sound of perturbations close to the 
source, we find, in contradiction to their original claims, that the presence of a su- 
perluminal mode cannot be avoided. This is a nontrivial result, which can only be 
seen by carefully taking into account the constraints that stability at large distances 
places on the theory, and the interplay between these conditions coming from in- 
finity and the action for perturbations near the source. In other words, this is not 
a local result which could have been derived from the knowledge of the behaviour 
near the source only. We have also showed that there will always be superluminal 
perturbations around a point source if only a Cubic Galileon is present. 

Our results physically arise from the link between the Vainshtein mechanism 
and superluminalities in typical Galileon theories. So long as one is considering 
theories that are ghost-free, with trivial asymptotic conditions at infinity and avoid 
quantum strong-coupling issues (fields with vanishing kinetic terms), these two ef- 
fects are intimately connected. One way to see this link is to note that the Vainshtein 
mechanism is inherently nonlinear, and so the behavior of the perturbations depends 
strongly on the source distribution present. Thus one expects to always be able to 
find backgrounds around which there are superluminalities. However, the connec- 
tion may be stronger than this: As we have shown, even in the case of a static point 
source with only a Quartic Galileon present, where the presence of superluminalities 
at large distances is not manifest, there are still inevitably superluminalities close 
to the source. 

We would like to emphasize however that the presence of superluminal modes 
is not enough to conclude that Galileon theories are inconsistent. As discussed in 
greater detail in Ref. [45], the Galileons still have their own causal structure. The 
crucial issue is instead whether or not closed time like curves can form. This would 
lead to violations of causality and the theory would be inconsistent. However, a 
Chronology Protection Conjecture for Galileon theories can be constructed, which 
states that it is impossible to form closed timelike curves without requiring energies 
that push the theory beyond its regime of validity. Further work could explore this 
conjecture in greater detail. 

We believe that the superluminalities are a crucial feature of Galileon theories. 
As shown in Ref. [38], the presence of superluminalities around some backgrounds 
is ultimately tied with the failure of the Galileon theory to have a Wilsonian com- 
pletion. It would be interesting to understand whether this aspect and the presence 
of a Vainshtein mechanism could however be tied to theories which allow for an 
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alternative to UV completion such as classicalization, [50, 51, 53]. 

We conclude by reviewing the only known way (so far) to have a Vainshtein 
mechanism and still avoid superluminalities. If the Galileon is not considered as a 
field in its own right, but rather as a component of anther fully fledge theory, one 
needs not to impose trivial asymptotic conditions at infinity. In massive gravity for 
instance, the Galileon field that appears in its decoupling limit is not a fundamental 
field. In such setups, it is then consistent to consider configurations for which the 
Galileon field does not vanish at infinity, so long as the metric is well defined at 
infinity (which does not necessarily imply Minkowski space-time). In such cases, 
we can thus have more freedom to fix the asymptotic boundary conditions for the 
Galileon field. A specific realization has recently been found 12 in [21], where the 
asymptotic behaviour is non trivial and the metric asymptotes to a cosmological one 
at large distances. These results do rely on the existence of a non-trivial coupling to 
matter of the form d^ird^T^ which naturally arises in Massive Gravity, [32]. When 
such non-trivial asymptotics conditions are considered, the results derived in this 
work are no longer valid and open the door for a way to find configurations which do 
exhibit the Vainshtein mechanism without necessarily propagating a superluminal 
mode around these configurations. Future work should consider the role of boundary 
conditions in the selection of viable configurations. 

Comments 

Part of this work was been derived as Paul de Fromont's Master thesis in the summer 
of 2011 and was presented in Ref. [52]. 

As a result to this work, the example and conclusions originally presented in 
Ref. [49] have been corrected, and a new arXiv version of [49] has been submitted, 
taking into account and summarizing the new analysis performed here. 

The results presented in this work are in agreement with that of Ref. [54]. 
Whilst the methods used in [54] are different, they reach the same basic conclusion 
that the combination of an active Vainshtein mechanism, a lack of ghost at infinity, 
and trivial boundary conditions at infinity inevitably lead to superluminalities in 
Multi-Galileon theories. 
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12 In a single Galileon realization. Notice that the presence of multiple Galileon fields would not 
be relevant there. 
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A Detailed analysis of the Special Case in the Quartic Galileon: 
Dominant First Order Corrections 

In this section of the appendix we will show that r — 5 > when Af > A|. 
We can write r — 5 as 



6b 01 m 

r — o = — — (1 — a) x 



1 / 3a — 1 2 

, — TTTj oioOoi - 2- -a 01 

Ajyi L8(3a — I) 2 \ 5a — 1 

1 



+ (1 - a)(5a - 1) C (5a - l)(3a - l)^ 01 



(A.l) 



We also use the notation A2 = aA 2 . Note that since ?/i,6oi > we must have 
3A 2 > A 2 . Thus we have 1/3 < a < 1, the upper bound comes from our assumption 
that Af > A|. 

Now the first term in the brackets of the expression r — 5 has the same sign as 

aoifroi - e«oi > ( A - 2 ) 

with < e < 1. This is positive because in order to avoid ghost instabilities at large 
distances from the source aio&oi — > and aio,&oi > 0, (see equation (5.1)). 
Meanwhile, the second term in the brackets is manifestly positive. Finally, the third 
term has the sign of — Cooi- 

So at this point our only hope of avoiding superluminalities is to consider a 
choice of parameters where — C«oi < 0. Now we will proceed to show that r — 5 > 
in this case as well. 

Note that in the limit a — > 1 with everything else fixed we have 

3 &oimC 2 

Also in the limit a — > 1/3 with everything else fixed we have 

2Afe, (l-ito)' >0 ' <A - 4) 



„^=^>0. (A.3) 



Now consider the function 



a(a) = 8(1 - a) (5a - l)(3a - 1 



,2 



«io^oi — 2 



(3a — 1J- 4 \ 5a — 1 



2 .2 1 



C - ~Tr 7T7^ rrCaoi 



(l-a)(5a-l)* (5a-l)(3a-l) 



(A.5) 



We can write this function in the shortened notation as tr(a) = cr + a\a + c^a 2 
with 

cx = 2a 2 ,! - ai 6oi + 8a iC + 16C 2 (A.6) 
a x = -8a 2 01 + 6a 10 6 i - 32a iC - 96C 2 (A. 7) 

a 2 = 6a 2 ! - 5a 10 6 i + 24a iC + 144( 2 . (A.8) 
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The sign of a(a>) is the same as the sign of r — 5 in the regime 1/3 < a < 1. Note 
that er(l/3), cr(l) > using the limits above. 

Note that 0" o ,o"i,er 2 do not have a definite sign, because a iC,C 2 > 0, but 
a oi — Oio^oi < 0. Therefore we need to investigate the behaviour of this function 
a (a) in more detail. 

Being a quadratic function a(a) has a single critical point (either corresponding 
to a maximum or a minimum) a cr i t . Given that <r(l/3), cr(l) > 0, we can avoid 
superluminalities if and only if 1/3 < a crit < 1 and simultaneously a(a crit ) < 0. 

Computing = yields for the critical point a cr i t 

a x Aa 2 m - 3a 10 b 01 + 16a iC + 48( 2 . . 

a C rit = -— = — • ( A - 9 ) 

2<7 2 0"2 

Plugging this back into the expression for a(a crit ) gives the following expression 

/ \ At l 2 \ a 01 + 8 a 0lC + 16C 2 — a 10^01 t\in\ 

v{a c nt) = 4(ai i - a i) ■ (A.10) 

0"2 

It is useful to consider 

1 - a crit = 2 < + 4ao 1 C + 48C 2 -a 1 oV 

0"2 

If a cr j t < 1 then this is positive. Similarly 

1 _ 1 6a§! + 24a iC ~ 4ai &oi c a i o\ 

®crit — 77 — 77 • l A - i/J J 

3 3 a"2 

If a cr u > 1/3 then this is positive 13 . 

We will now show that we cannot simultaneously satisfy all the criteria that we 
need to satisfy to avoid superluminalities. We consider four cases which will exhaust 
all possibilities: 

Case 1: a 2 = 
In this case we have 

a(a) = (To + o'lOt 

Since <r(l/3), cr(l) > we know that a{a) > in the whole interval 1/3 < a < 1. 
Case 2: a 2 < 

Consider a(a cr it). If we assume that a 2 is negative, then we can avoid superlumi- 
nalities if and only if the numerator of u{a cr it) is positive. 

However the condition that cr 2 is negative means that aoi&oi > f^oi + y a oiC + 
C 2 , which implies that 

num of a{a crit ) = a 2 01 + 8a iC + 16<C 2 - aio&oi < -^(aoi - §C) 2 • (A.13) 
So we cannot avoid superluminalities in this case either. 



144 
5 



13 Whcn we write num of o~(a cr it), we mean + 8aoi£ + 16£ 2 — aio^oi by that, ie we are 
ignoring the uninteresting factor of 4(aio&oi — a oi 

) > 0. 
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Case 3: a 2 > 0, C > 

Again we consider a(a crit ). We now assume that cr 2 is positive, so we need to 
check if numerator of a(a cr i t ) can be negative if we also assume that a cr n < 1, ie 
aio^oi < «oi + 4aoiC + 48£ 2 , and also that a cr u > 1/3, ie aoi&oi < f^oi + SooiC 
The inequality a cr u < 1 tells us that 

num of a(a crit ) = + 8a i( + 16( 2 - a 10 b i > 4a iC-32C 2 = 4C(a i-8C) (A. 14) 

and the inequality a crit > 1/3 tells us that 

num of a(a crit ) = a 2 01 + 8a 01 ( + 16( 2 - a w b i > -7^01 + 2a oiC + 16( 2 (A. 15) 

Now let's take ( > 0. The first inequality then implies we need aoi — 8( < to avoid 
superluminalities. So we set a i = 8(e for < e < 1 (if e < then — a iC > 0). 
Then the second inequality becomes 

num > 16C 2 (1 + e - 2e 2 ) = 16C 2 (1 - e)(l + 2e) > (A.16) 

So also in this case we are forced to have superluminalities. 

Case 4: a 2 > 0, ( < 

Now we take ( < 0. The first inequality then implies we need aoi — 8^ > to avoid 
superluminalties. However note that both aoi and £ are negative. 

So we set aoi = 8(e for < e < 1. Then the argument is exactly the same as 
above, and that concludes our set of possibilities. In conclusion there is no possible 
way to avoid superluminalities near the origin, even if one had been so lucky as to live 
in a theory with specifically tuned coefficients for which the first order corrections 
near the origin vanished. Our result is thus generic: superluminalities are always 
present near the origin if the field is to be trivial at infinity and stable both at small 
and large distances. 
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